Localization and shock waves in curved manifolds for the Gross-Pitaevskii equation 
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We investigate the dynamics of a Bose-Einstein condensate in a progressively bended three dimen- 
sional cigar shaped potential. The interplay between geometry and nonlinearity is considered. At 
high curvature, topological localization occurs and becomes frustrated by the generation of curved 
dispersive shock- waves when the strength of nonlinearity is increased. The analysis is supported by 
four-dimensional parallel simulations. 



The effect of geometry on wave propagation has been 
considered since early developments of the theory of 
sound, for example, in the vibrations of membranes with 
different shapes [1]. Geometry largely affects energy ve- 
locity, wave transformations upon propagation, and vari- 
ous classical and quantum undulatory phenomena, either 
in constrained geometries, or at an astrophysical scale; 
renowned examples are the Einstein lensing effect (H , the 
Hawking radiation [3,] ,01" the Unruh effect [j] 

Albeit this old history, the link between geometry and 
waves is continuously fascinating many researchers, and 
is currently driving_^ import ant applications as transfor- 
mation optics jsl, lil, curved and twisted waveguides, 
[5I, and analog of gravity in Bose-Einstein condensa- 
tion (BEG) [8, 9] or nonlinear optics [10-12]. It is also 
known that nonlinear waves may be affected by geomet- 
rical constraints, resulting in interesting phenomena, as 
solitons in curved manifolds 0, ll^ . 

A striking effect of curvature is the topologically in- 
duced localization, originally considered in [15]: when the 
wave is constrained on extremely deformed surfaces, or 
lines, propagation may be inhibited and energy is trapped 
in the regions with the highest curvature. The way this 
kind of geometrical localization competes with nonlin- 
earity, and the kind of dynamic effects resulting from a 
nonlinear response sufficiently strong to overcome topo- 
logical bounds is largely unconsidered. 

If geometrical localization occurs in a reduced dimen- 
sionality (e.g., a curved surface in a three dimensional 
space), when nonlinearity is very effective the whole 
three-dimensional (3D) space becomes involved, and any 
treatment based on nonlinear wave equations with re- 
duced dimensionality may be questioned. This is a key 
difficulty in this problem; any theoretical prediction must 
be numerically tested by using 3D simulations. 

In addition, for high nonlinearity, shock waves (SW) 
are generated [16] . In recent years^W have been largely 
considered in optics and BEG |17H2l| ; SW originate from 
singular solutions of the hydrodynamic reduction of the 
nonhnear Schroedinger equation (NLS) regularized by os- 
cillating wave fronts, named undular bores, or dispersive 
SW (D-SW). To the best of our knowledge, the effect of 
a curved space on DSW is un-explored. 

In this Letter, we investigate the way a geometrical 
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FIG. 1. (Color online) (a) Potential profile for the piecewise 
case (dashed, R — 0.5) and for the parabolic cases {y — kx^ 
with R = 1/2A; = 1,2,3); (b) geometrical potential Vg for 
the cases in panel (a); (c) shape of the bound state for the 
parabolic potential for various powers (PR goes from 10~^ 
to 2.6); (d) nonlinear eigenvalue versus power (left axis), and 
localization length (right axis), wi R^ = 6.25 X 10"^ 



localization is frustrated by a defocusing nonlinearity, 
we study this effect by one-dimensional (ID) theoretical 
analysis and 3D+1 simulations of the Gross-Pitaevskii 
(GP) equation, and report on DSW in curved potentials. 

The model — The adimensionalized GP equation is 

i^t = - V'^ + V^r)^ - xl^l'^. (1) 

In ([1]), X = determines the sign of the nonlinearity, 
and N = J |^p(i^r. We consider a curved potential, 
sketched in figure [U Vsd = Vi{qi) + ^±(^2,^3) where 
qi = q is the longitudinal coordinate along the arc, and 
q2^3 are the tran5i;er5e coordinates. [15] V± determines the 
transverse confinement along an arbitrary curve, whose 
curvilinear coordinate is given by q. Vi{qi) is a weak lon- 
gitudinal trapping potential, whose effect becomes neg- 
ligible for strong curvature, as detailed below. Specif- 
ically, we have a parabolic potential V_\_ = w{q2 + ^3), 
"^1(^1) = '^1^1 with w » wi. Being = w/wi, the ID 
reduction holds true as ^ 00 [15], letting 



^(<7i,^2,g3) = ((72,^3) V^(<7i)exp(-z£;xt) (2) 
'^{Q2^q3)dq2dq3 = 1, and Et the trasverse part 



with J 

of the eigenvalue. The transverse localization length is 



It = / lV^±(<72,g3) 
dition J\i^{q)\'^dq 



dq2dqs^ and the normalization con- 
= 7V/J^ = P. For the parabolic 
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FIG. 2. (Color online) Snapshots of the two dimensional section of the wave- function in the bending plane xy for different 
curvatures in the linear case (x = 0? ^ — 1-5); (a) k — 0.5, (b) k — 1 and {c) k — A 



potential above, one has for the linear ground state 
= v'27^exp(-^|//2 - ql/ll) with It = (4/^)^/^ 
and Et = 2y/w. The resulting ID GP equation with 
geometrical potential Vq is given by 



(3) 



where V = Vi -\-Vg^ and Vg is expressed in terms of the 
local curvature K{q): 



Vaiq) 



K\q) 



(4) 



The bound state can be determined by the normalized 
one-dimensional equation 



(5) 



Linear geometrical localization — For x = and Vi = 0, 
we consider a step- wise curvature [l5j , as in figure [TJi, 
dashed line: Vciq) = for \q\ > R and Vciq) = 
— l/(4i?^) for \q\ < R with R the minimum radius of 
curvature. This piecewise potential allows simple ana- 
lytical results: there is a single bound state, which is 
exponentially localized with respect to q and exponent 

with E -0.1/i?2 < 0. 
For a bended cigar shaped trap (FigUJi), with profile y = 
kx'^, we have 



Vaiq) 



K{q? 



[l + AkH{qY 



(6) 



where x{q) is given by the inverse of Akq = 
2A:xvTT4Fx2 + sinh"^(2A:x). The maximum curva- 
ture is 2k = 1/R. As for the stepwise potential, a 



unique eigenstate with negative energy exists with E = 
— as numerically calculated. An analytic approx- 
imation is found by replacing Vg by a Dirac 5 with the 
same area, this furnishes E = — l/(9i?^) and decaying 
tail (X exp(— Ig^l/Zc), with localization length Iq = 4i? of 
the order of R. 

When also considering the presence of a longitudinal 
trapping potential Vi{q) 

sian bound state (p = 

h = V2/i 



wiq^ 



the associated Gaus- 



y V2/7rP//iexp(-^V/^) with 

/w^^^^ may mask the exponential localization 
as far as /i is comparable with the local radius of curva- 
ture; geometric localization appears for Iq « h- 
Nonlinearity and curvature — For the piecewise curva- 
ture, in the presence of defocusing nonlinearity (x = —1), 
the Thomas-Fermi approximation allows to write cp'^ = 
E - l/(4i?2) for \q\ < ttR/2 and = elsewhere 0. 
The normalization J cp'^dx = P provides the power de- 
pendent nonlinear eigenvalue E = — l/(4i?^) + P/(7ri?), 
which shows that, for P = Pth = 7r/(4i?), E changes 
sign. For P > Pth^ E > 0, and the wave is delocalized: 
the nonlinearity destroys the geometrical localization. 

For the parabohcally bended cigar-shaped potential for 
X = —1, we numerically solve the nonlinear bound-state 
equation ([5]) by a pseudo-spectral Newton-Raphson ap- 
proach, also including the longitudinal potential Vi = 
wiq'^. In figure [it we show the calculated profiles for 
different powers; the results can be expressed in terms 
of variables suitably scaled by R without loss of general- 
ity. When increasing P, there is a transition from an ex- 
ponentially localized state to the Gaussian ground state 
of the longitudinal potential; correspondingly E changes 
sign, from negative to positive (Fig(T]i), and the local- 



3 




nonlinearity N 



evolution t 



0.5 1 
evolution t 



FIG. 3. (Color online) (a) Linear case (x = 0), localiza- 
tion length versus time for various values of the curvature k 
{k — 0,0.5,1,2,4); (b) nonlinear case {k — 1^ x — "1)? lo" 
calization length versus time t for various number of atoms 
{N = 10,20,30,100,500). The localization length is calcu- 
lated as the three-dimensional inverse participation ratio. 



ization length increases. When Vi = a transition to a 
delocalized state is found but not reported. Note that 
the power for delocalization scales as the inverse of the 
minimum radius of curvature, i.e., the more the curva- 
ture, the higher the strength of nonlinearity needed to 
delocalize. This in agreement with the simple theoreti- 
cal calculations for the stepwise potential reported above. 
We numerically find Pth = 0.5/ R k. 
Numerical results in the linear case — We solve the 
3D+1 GP equation by a pseudo-spectral parallel algo- 
rithm. We first consider the case without nonlinearity, 
and show the way the geometrically induced localiza- 
tion appears in the spreading of an initially localized 
wave packet in a cigar shaped potential when increas- 
ing curvature. We consider an initial wavefunction that 
matches the transverse size of the linear ground state for 
the transverse potential V± = wq^^ with w = 100; the 
initial condition is ^ = exp(— g^//^) exp(— — Qs/It) 
with li = It = (4/k;)^/^ = 0.44, other parameters are 
= 10^ and Vi = wiqf and Wi = 10~^, the initial size 
is indicated by the circle in figl2^. 

When starting from low bending {k < 1), as the ini- 
tial longitudinal length is much smaller (/i = 0.44 in our 
units) than the size of the fundamental state (/i = 8), 
the wave spreads in the q direction, as shown in figure 
[2^. When the curvature increases, the wave displays a de- 
gree of localization, corresponding to the excitation of the 
topologically trapped state, as shown in Figl2l3,c. The 
related localization length is of the order of the radius 
of curvature. FiglSt^ shows the localization length cal- 
culated by the inverse participation ratio for an increas- 
ing curvature; at high curvature the localization length 
asymptotically becomes nearly time-indepedent. 
Hydrodynamic approximation — In the presence of a de- 
focusing nonlinearity, the geometrical localization is op- 



posed to the tendency to spread. Following the analysis 
above, one expects a power threshold for the nonlinear 
frustration of the topological trapping; we show in the 
following that this is mediated by the generation of DSW. 

We consider an initial wave-packet with Gaussian 
shape and waist of the order of the curvature, as in Figl2l 
We investigate the expansion of this condensate in the 
presence of a defocusing nonlinearity and of the geomet- 
rical potential. The theoretical analysis can be developed 
in the hydrodynamical approximation j23[. Within the 
ID model, we consider the evolution of an Gaussian beam 
ip{t = 0) = exp(— 7^^/2) with P = v^/e. Introducing 
T = t/e, T] = eg/>/2, the scaled equation reads as 



ietlJr 



(7) 



In the limit e ^ 0, the occurrence of an hydrodynamic 
shock is signaled by collapsing trajectories r]{r) solutions 
of the equation |2^ 



df]'^ _ dVhydro _ d 



dr^ 



df] 



df] 



(8) 



The dynamics is determined by the potential Vhydro — 
VG-\-Vp, with Vp = exp(— 7^^) given by the initial Gaussian 
wave-profile. In the absence of curvature Vg = and 
the shock is determined by an inverted well Gaussian 
potential Vp] in this case it is well known that the wave 
develops a couple of symmetrical lateral shocks [24^. In 
the presence of curvature, as < and Vp > 0^ the 
geometry seems to prevent the shock, however, due to 
the specific properties of the topological potential, this 
effect is negligible and curved DSW are expected. 

Indeed, following Eq.(|6]), the peak value of the geo- 
metric potential is k'^ and its width in the scaled variable 
T] is Wr^ e//c, while the width of Vp is of the order of 
the unity. There are hence two possibilities: (i) when 
starting from a wavefunction with size larger than the 
radius of curvature, the width of the potential well w^^ 
is much narrower than the initial wave profile, and Vg 
does not prevent the occurrence of the lateral shocks; (ii) 
when considering an initial condition with size compara- 
ble with the radius of curvature (k;^ = 1, e = /c), the 
maximum value of Vg, k'^, is much smaller than the uni- 
tary peak value of Vp and the shock occurs when e << 1, 
i.e., when P » ^JH. 

The geometrical potential and the related topological 
localization hence do not prevent the formation of DSW, 
which occur along the curved manifold. One may argue if 
the approximations underlying the reduced model Eq.(|3]) 
are still in the presence of strong nonlinearity, this is 
confirmed in the simulations reported in the following. 
Dispersive shock waves in curved manifolds — Here we 
consider the evolution of an initially localized wavefunc- 
tion with nonlinearity; we study the 3 + ID evolution 
when increasing power. A representative series of snap- 
shots are given in figure HI we show sections in the 
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transversal plane xy where curvature occurs, the wave 
stays confined by the trap in the xz plane. The occur- 
rence of the characteristic oscillation associated with the 
generation of the undular bore in the strong nonlinear 
regime is evident. The shock formation is due an over- 
whelming nonlinearity that masks the geometrical trap- 
ping. The trend of the localization length versus time for 
increasing power is in figure [SJd, which also shows oscil- 
lations typical of the undular bore. 

Conclusions — We investigated the interplay between 
geometrically induced wave-localization and nonlinearity 
for the Gross-Pitaevskii equation. For a defocusing non- 
linearity, we determined the power needed for the tran- 
sition to a delocalized regime in terms of the curvature. 
We have shown that at the frustration of the geomet- 
rical trapping, the nonlinearity induces hydrodynamical 
phenomena. The theoretical analysis is supported by nu- 
merical simulations in four dimensions. Our results show 
the way complex manifolds in Bose-Einstein condensa- 
tion, and in related fields as nonlinear optics, may systain 
novel dynamical effects, as curved dispersive shock waves, 
and open the road to a variety of related analyses on the 
competition between geometry and wave-localization. 
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